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The method used at the National Bu- 
reau of Standards for the calibration of 
strontium-90 + yttrium-90 beta-particle 
ophthalmic applicators in terms of ab- 
sorbed dose to water, is described. The 
method involves measurement of ioniza- 
tion density at the applicator surface 
with an extrapolation chamber, correc- 
tion for the difference in backscatter be- 
tween the collection electrode and 
water, and application of the Bragg- 



Gray equation. The calibration obtained 
is an average over the active surface of 
the applicator. The overall uncertainty 
of the surface calibration is about ±15 
percent. 

Key words: absorbed dose; calibration 
service; electron dosimetry; ophthalmic 
applicator; surface dose. 

Accepted: August 27, 1986 



1. Introduction 

This paper describes the calibration of strontium- 
90 4- yttrium-90 beta-particle ophthalmic applica- 
tors in terms of absorbed-dose rate to water at the 
surface of the applicator. The calibration is made 
with an extrapolation ionization chamber, an in- 
strument invented by Failla [l] 1 which has since 
been shown by several investigators [2,3] to be par- 
ticularly well suited to the determination of ab- 
sorbed-dose rate for beta-particle sources. The 
techniques of measurement have been sufficiently 
refined to permit the successful international com- 
parison of several sources at distances of 20 and 30 
cm, to within one percent [4]. 

For ophthalmic applicators, the quantity of inter- 
est is absorbed-dose rate at the applicator surface. 

About the Author: J. S. Pruitt is with NBS' Ioniz- 
ing Radiation Division, Center for Radiation Re- 
search. 

1 Figures in brackets indicate literature references. 



The problems introduced by reducing the source 
distance from 20 cm to zero make calibration un- 
certainties of about 15 percent realistic. Within this 
limit, it has been found that a reliable surface cali- 
bration of an ophthalmic applicator can be ob- 
tained with an extrapolation chamber. The 
extrapolation chamber technique has been used at 
the National Bureau of Standards (NBS) since 1977 
to calibrate about 10 ophthalmic applicators per 
year. Alternative methods, using thermolumines- 
cent dosimeters (TLD) have also been described 
[5]. The use of TLD 's offers a simpler and less elab- 
orate approach, but it requires calibration of the 
TLD material in another radiation source, and was 
not considered an appropriate method for use at 
NBS. 

The general plan of the calibration is to measure 
ionization current per unit mass of air with the ap- 
plicator at several distances from the center of the 
extrapolation chamber air gap, and to extrapolate 
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to zero distance to obtain a value at the surface of 
the applicator. The following three sections are 
concerned with this determination of extrapolated 
current per unit mass of air. This is followed by a 
section describing the calculation of the absorbed- 
dose rate, and a final section about the errors inher- 
ent in this calibration method. 



2. The Extrapolation Chamber 

The extrapolation chamber is essentially the 
same as that described by Loevinger and Trott [6], 
with several collecting electrodes which can be in- 
terchanged and a thin aluminized-plastic-film high- 
voltage electrode. The electrode spacing can be 
changed by rotating the high-voltage-electrode 
support, at the rate of 0.635 mm per revolution. 
The relative position of the electrodes and an oph- 
thalmic applicator under study are shown in figure 
1 . The particular collecting electrode shown has a 
diameter of 28 mm and is surrounded by a guard 
electrode with an outside diameter of 50 mm. Both 
guard and collecting electrodes are made of 8 mm 
thick A- 150 conducting plastic, electrically isolated 
with 30-fxm plastic film held in place with epoxy 
cement. The high-voltage electrode is aluminized 
plastic with a diameter of 80 mm and a thickness of 
0,7 mg/cm 2 . 

The collecting electrode shown in figure 1 has a 
much larger area than the active area of all ordi- 
nary ophthalmic applicators. This type of collect- 
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Figure 1-Schematic diagram of ophthalmic applicator and cross 
section of extrapolation chamber electrodes. 



ing electrode is used to determine the ratio (I/d) , 
where I is the ionization current generated in the 
air volume defined by the collecting area, A, and 
the electrode separation, d. The ratio I/d is a func- 
tion of both d and the applicator distance, z, and 
the zero subscript refers to the limiting case where 
both d and z are zero. This limiting value of I/d is 
a characteristic of the source and of the electrode 
material and is independent of A, if A is sufficiently 
large. 

A probe electrode of smaller area is used to in- 
vestigate the variations in dose rate across the sur- 
face of applicators. This probe is similar to the unit 
shown in figure 1, but the diameter of the collect- 
ing electrode is only 0.6 mm, and both collector 
and guard are made of graphite. A 0.4 mm thick 
ring of epoxy between collector and guard results 
in an effective probe diameter of about 1 mm. With 
this probe, both the local activity variation and the 
effective source area are investigated by changing 
the x and y coordinates of the applicator (fig. 1), 
for small but fixed values of z and d. The applicator 
stand rests on a microscope stage to permit accu- 
rate variation of both x and y. Accurate adjustment 
of z is made with a barrel micrometer. 

Extrapolation chamber currents include contri- 
butions from (1) ions generated in the collecting 
volume of the chamber, (2) electrons coming to 
rest in the collecting electrode, and (3) back- 
ground. The first of these is the desired quantity. It 
is separated from (2) and (3) by measuring both I + 
with positive high voltage (a positive number) and 
I - with negative high voltage (a negative number). 
Since only (1) changes sign with the high voltage, 
the ion current is (I + — 1~)/2, a quantity which is 
proportional to the electrode separation, to a first 
approximation. The quantity (I + + I _ )/2, called the 
"parasitic" current by Hillion et al„ [7] includes 
contributions (2) and (3) and is almost independent 
of electrode separation for a given applicator dis- 
tance. 

Figure 2 shows the variation of both the ion and 
parasitic currents across the face of the most asym- 
metric applicator that has been studied in this labo- 
ratory. The measurements were made with the 
probe electrode, using an electrode spacing of 0.05 
mm and an applicator distance of 0.4 mm. 2 The 
asymmetry of the ion current is presumably caused 
by an uneven distribution of activity across the ap- 
plicator face. The tails of the ion current distribu- 
tion are caused by beta rays that leave the 
applicator with large velocity components perpen- 
dicular to the beam direction. The tails rapidly be- 

2 In measurements made with the 28 mm collecting electrode, 
the parasitic current is only 5-10 percent of the ion current. 
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Figure 2- Variation of ion current and parasitic current across 
applicator SN 0136. 

come more prominent as the applicator distance z 
increases. 

The graph of parasitic current in figure 2 shows 
evidence both of electrons coming to rest in the 
collecting electrode (the central peak of the distri- 
bution), and background (the negative tails). The 
central peak of the parasitic current is more sym- 
metrical than the ion current peak, probably be- 
cause of multiple scattering of the primary beam in 
the collecting and guard electrodes, which tends to 
round off the edges of the distribution. The back- 
ground is a leakage current of about 10-50 fA. It is 
unaffected by high-voltage polarity changes, as 
stated earlier. The parasitic currents are generally 
ignored. 



3. Determination of Source Area 

Initial inspection of an applicator is made with 
the probe electrode, using as small an electrode 
separation and applicator distance as can be toler- 
ated. The applicator face is mapped at points sepa- 
rated by 2 mm, over a grid of sufficient extent for 
the ion current to drop below half of its maximum 
value along all tracks. This crude map permits cal- 
culation of the source area and provides a qualita- 
tive picture of the source symmetry and activity 
distribution. On occasion, when more detail is re- 
quired or when the source diameter is found to be 
only a few millimeters, a 1 mm or even a 0.5 mm 
grid spacing is used. 



The source area is taken to be the area within the 
locus of half-maximum-current points on such a 
map. It is determined from a plot of the measured 
current as a function of radius from the center of 
the distribution, as shown in figure 3 for a particu- 
lar applicator, For most applicators, the current de- 
creases monotonically with increasing radius. 3 A 50 
percent radius, R, can be read directly from the 
graph, and the area calculated from this as 7rR 2 . 
This is a close approximation to the real area, ttR 1 , 
if the chamber spacing, d, and source distance, z, 
are both sufficiently small. A rough estimate of un- 
certainty for a given d and z can be obtained from 
the following equation describing the ionization 
current per unit volume as a function of position: 



dI V 1 
^ = K.log e 



3 For some applicators, the probe electrode map shows re- 
gions of high activity near the edge of the distribution. These 
hot spots make the determination of area very uncertain and at 
the same time bring the medical use into question. We decline to 
calibrate such sources and recommend against their use as oph- 
thalmic applicators. 




RADIUS, mm 

Figure 3-Currents measured with the probe electrode across the 
face of SN D-402, plotted as a function of radius, measured 
from the center of the current distribution. The curve was 
drawn freehand. 
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(1) 



where in terms of th e coor dinates of Figure 1, 
fi = r/R, tj = z/R, r=Vx 2 -t-y 2 , and K is a constant. 

This equation was derived assuming that the ac- 
tivity is uniformly spread over a circle of radius R 
on the source surface, that the beta-particles are 
emitted isotropically from the source surface, that 
their paths are straight, and that energy deposition 
is uniform along these paths. The derivation is 
based on a 1/r 2 distance variation and consequently 
current per unit volume becomes infinite as the dis- 
tance approaches zero. If the current/volume is 
normalized to its value on the central axis of the 
source (jx=Q), the normalized current is well be- 
haved for all distances (that is, for all values of tj). 

Figure 4 is a plot of the normalized currents for 
several values of T). For each curve, the apparent 
radius, R, is the product of the true radius, R, and 
the value of (i where the curve crosses the half- 
value line. For tj=1, z=R, and 7rR 2 exceeds 7rR 2 
by a factor of about 2. For 17=0.2 (not shown) irR 2 
exceeds 7rR 2 by about 4 percent, and for 17 =0.1, 
the discrepancy is less than 1 percent. In a typical 
case, d=0.2 mm, z = 0.12 mm, R = 4 mm, and 
17=0.03, for which 7rR 2 /VR 2 = 1.001. 

Figure 5 is a plot of probe ion current vs. radius 
for two applicators and matching source profiles 
calculated from eq (1). The parameter z was as- 
signed the approximate value 0. 1 mm in the calcu- 
lation, and values of R=3.18 and 4.35 mm were 
read from the graph. All four curves have been 
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Figure 5-Measured probe profiles of applicators SN 0258 and 
SN 157, in comparison with theoretical predictions for 
uniformly distributed sources. 

normalized at the origin. Equation (1) is not an ac- 
curate description of either source profile, but the 
comparison with SN 0258 is consistent with the 
assumption that the activity in this source is spread 
uniformly over an area of 4.35-mm radius. The 
small differences between theory and measurement 
appear to be easily explained. The differences near 
r=R can be attributed to the finite size of the probe 
electrode. The differences at larger radii probably 
arise from the fact that the activity in SN 0258 lies 
below the applicator face, so that the elemental 
sources cannot be treated as isotropic at large an- 
gles. This will cause the real distribution to fall off 
more rapidly at large radii than the theoretical dis- 
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Figure 4-Theoretical probe profiles 
for five different source distances, 
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tribution. The activity in SN 157, in contrast, 
pears to be highly peaked at the center.* 



ap- measurement. Since the currents have been cor- 
rected to the reference temperature and pressure, 
the relevant density is 1.197 kg/m\ 



4. Determination of Average Current 
Per Unit Mass of Air 

After mapping the source and measuring its area, 
the probe electrode must be replaced by the larger 
collecting electrode to determine the current per 
unit electrode spacing (I/d)t>. This is determined by 
measuring current, I, for five different air-gaps, 
d = 0.5, 1, 1.5, 2, 2.5 mm, with the source distance 
z=d/2 in each case. These measured currents are 
corrected (1) to a reference air temperature and 
pressure, 22 °C and 101.3 kPa; (2) for ion recombi- 
nation and diffusion, using the formulas and con- 
stants summarized by Bohm [8J; and (3) errors in 
the extrapolation technique, discussed below. 5 The 
corrected current is a non-linear function of d, be- 
cause of poor geometry, and extrapolation to zero 
d is imprecise. The corrected current is fitted to a 
cubic equation in d constrained to pass through the 
origin. The ratio 1/d is then a quadratic function 
that can be written: 



I 

d 



@. 



(t + ad-rtfd 2 ) 



(2) 



-0.1 mm -1 and £=0,01 



In a typical case, a% 
mm -2 . 

The (I/d)ti ratio from eq (2) depends on the ex- 
trapolation path to a minor extent. For example, if 
the ratio z/d had been maintained at 1 rather than 
at 1/2, the measured (I/d) ratio would have been 
slightly smaller. Ideally, the ratio should be mea- 
sured varying only d, with z—0 throughout. The 
closest one can come to this ideal situation is to 
extrapolate I/d to zero d for several fixed values of 
z and then extrapolate the zero-d limits to zero z 
independently, a procedure which increases the 
calibration time by a large factor. For applicator 
SN 0258, it was found that the ideal (I/d) is 2 per- 
cent smaller than the value obtained with z/d=l/ 
2, This is the origin of the correction factor (3) 
above, which is set equal to 0.98. 

For a given source, the average current per unit 
mass of air is the ratio of (I/d) to source area, A„ 
divided by the air density, p, at the time of current 



4 A closer theoretical fit may be obtained by using ij = 0.5 for 
SN 157, but this, requires the assumption of an unrealistically 
large z. 

5 No corrections were made for the fail deformation effects 
discussed by Loevinger. These effects are too small to be de- 
tected with the Mylar foils used. 



5. Calculation of Absorbed Dose 

The average surface absorbed-dose rate in water 
is given by a variant of the Bragg-Gray equation: 






back "Vfc-Ll 



(3) 



where W/e is the average energy_ required to pro- 
duce unit charge in ambient air, (S/p)™ ,er is the ra- 
tio of the average collision mass stopping power of 
water to that of air, k^v. corrects for a difference in 
backscatter between the collecting electrode and 
water, and k M corrects for attenuation or buildup 
in the high-voltage electrode. 

The value, W/e = 33.7 J/C, is used in eq (3) [9], 
The stopping power ratio was calculated with un- 
restricted stopping powers [10], using two different 
source spectra. The first spectrum was based on 
unfiltered theoretical spectra for strontium-90 and 
yttrium-90 [11], modified by the inherent filtration 
in a particular source, SN 0258, and yielded the 
value {S/p)£r tE! = 1.123. The second spectrum was 
measured with a lithium-drifted germanium detec- 
tor for a strontium-yttrium source filtered by 50 
mg/cm 2 of silver and 0.1 mm of steel and yielded 
the value (S/pJ^ CSI = 1.124, the value used in our 
calibrations. 

The backscatter factor was estimated from pub- 
lished formulas for backscatter probability as a 
function of atomic number and beta-particle energy 
[12]. Using the chemical formulation for the col- 
lecting electrode material, and the second source 
spectrum described above, the value k ba - k = 1.010 
was obtained. 

The foil attenuation correction factor, 
kfcj| = 1.003, was determined from measurements 
with extra thicknesses of foil. 

In summary, the four right-hand factors in eq (3) 
may be combined to form a single factor, 38.4 
Gy-kg/C. Multiplication of this factor by the ratio 
of current per unit mass determined with the tech- 
niques of the previous section gives average ab- 
sorbed-dose rate to water over the active surface of 
the applicator. 



°The value (S/p)S? lcr =1.15, based on earlier stopping-power 
tabulations, had been used until recently. 
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6. Calibration Uncertainty 

There is some uncertainty associated with each 
term of eq (3). Those determined by extrapolation- 
chamber measurements with the 28 mm diameter 
collecting electrode, such as (I/d) and k foi |, have 
statistical uncertainties of the order of 0. 1 percent. 
These are much smaller than the estimated uncer- 
tainties and are usually ignored. The estimated un- 
certainties are listed in table 1. Each is intended to 
be comparable to a standard deviation, in that the 
chance that the true value of the quantity lies 
within plus or minus the stated uncertainty is in- 
tended to be about 2 out of 3. 

The instrumental uncertainties include voltmeter 
readings, capacitance measurements, and pressure 
and temperature determination. The uncertainty in 
W/e comes directly from reference 9. The stop- 
ping-power-ratio uncertainty includes uncertainties 
in the stopping-powers themselves plus uncertain- 
ties in the beta-particle spectra. The uncertainty in 
(I/d) is an estimate of the reliability of a cubic fit 
to an extrapolation curve with z = d/2, and the cor- 
rections (1), (2), and (3) in the text. The uncertainty 
in surface area is a standard deviation, based on 
repeat measurements for several sources, using the 
normal 2-mm grid. 

When the uncertainties are combined in quadra- 
ture to obtain a combined uncertainty, the result is 
7 percent. This number is intended to have the ap- 

Table 1. Estimated uncertainties. 



proximate character of a standard deviation. Dou- 
bling it gives an overall uncertainty of 14 percent 
which can be interpreted as having the approxi- 
mate character of a 95 percent confidence limit. 
For convenience, this uncertainty is rounded off to 
15 percent. 

The surface area uncertainty is dominant in table 
i and is the only one which can be appreciably 
reduced at present. It can be reduced to about I 
percent by mapping a normal source with a 1-mm 
grid, but this step is so much more time-consuming 
that it is not offered as an NBS calibration service. 
In this case, however, the combined overall uncer- 
tainty would be about 4 percent. This would lead 
to a stated uncertainty of 8 percent with the ap- 
proximate character of a 95 percent confidence 
limit. 

As a more practical alternative, it has been sug- 
gested that the source surface area could be deter- 
mined with greater accuracy with radiochromic 
dye films. 
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